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INTRODUCTION 

This section presents techniques for the optimal selection and 
location of instruments in process grassroot designs as well 
as in retrofit scenarios to perform data gathering for control 
and monitoring, as well as fault detection and online optimi- 
zation. The main goal is usually to minimize the instrumen- 
tation investment cost subject to performance constraints. 
Some of these constraints are as follows: 

1. The network and the associated software should be 
able to obtain estimators of measured and unmeasured 
variables with certain accuracy (given precision and 
bias free) for monitoring and control purposes. 

2. Malfunctioning instrumentation should be able to be 
detected and its readings filtered through an algorithm 
that makes use of all measurements. 

3. Unsafe and/or faulty operating conditions and their 
origin should be able to be detected. 

4. Conditions and events that can impact efficiency and 
quality of products should be identified. 

In recent years, a new paradigm is emerging. In this 
paradigm, it is recognized that data and the information 
obtained from data processing have a certain economic 
value. Thus, optimal networks are those that maximize the 
difference between the aforementioned economic value 
and the cost. 

When data acquisition is solely supported by instrument 
measurements, each particular value of a variable is directly 
associated to the instrument that measures it. Very precise 
and failure-proof instrumentation, albeit also costly, would 
suffice to obtain precise and reliable estimates. However, 
since the probability of a sensor failure is finite, redundancy 
is usually used as a means of guaranteeing data availabil- 
ity. With redundancy comes discrepancy among measure- 
ments. If the noise associated with the signal is small 
enough, these discrepancies can be ignored, but normally 
these readings have to be reconciled. Reconciliation also 
brings improved precision of the estimates. In turn, redun- 
dancy is classified as: 

• Hardware Redundancy : when two or more sensors are 
used to measure the same variable. 


• Software or Analytical Redundancy : when a set of 

values of different variables is supposed to satisfy a 

mathematical model. 

A simple example of analytical redundancy is a unit with sev- 
eral input and output streams. As an example in this article, if 
one measures the flow rate of all streams with a single instru- 
ment per flow measured, there is no hardware redundancy. 
However, since the sum of the input flow rates has to be equal 
to the sum of the outputs, redundancy is analytical. In other 
words, there are two estimates available for each flow rate, 
one from its direct measurement and the other obtained using 
the material balance equations, which constitutes a conflict 
that needs to be resolved. 

Data reconciliation and gross error detection are tech- 
niques that help obtain precise estimates and at the same time 
identify instrument malfunction (tasks (a) and (b)). While 
data reconciliation relies on analytical constraints (typically 
material and energy balances as well as other thermody- 
namic relations) and mostly on least square estimation, other 
techniques used to obtain estimates, like Kalman filtering, 
pattern recognition, neural networks, principal component 
analysis, partial least squares (PLS), among others, rely 
on different concepts. There are several books that outline 
many aspects of data reconciliation [12,53-55,67,82,95], In 
addition, hundreds of articles have been devoted to the prob- 
lem. Kramer and Mah [43] discussed how good estimates 
of data can be obtained by means of models and it is for 
this reason that these activities have been named model- 
based monitoring. Thus, the field of sensor network design 
and upgrading has traditionally relied on model-based con- 
cepts simply because it is the reverse engineering problem 
of monitoring. 

Instrument malfunction is a term that relates to situa- 
tions that range from miscalibration or bias to total failure. 
In the absence of redundancy, miscalibration or bias cannot 
be detected (unless the deviation is so large that it becomes 
obvious). Redundancy, and especially software or analyti- 
cal redundancy, is the only way to contrast data and deter- 
mine possible malfunctions of this sort. This approach was 
named gross error detection by one group of researchers 
[18,32,33,56,81,83] and sensor/data validation or signal 
reconstruction by another group [34,78,89,90] and many 
others. Data reconciliation and gross error detection are 
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also illustrated in detail in section Data Reconciliation and 
Software Methods for Bias Detection. 


MEASURED AND KEY VARIABLES 

Traditionally, an a priori selection of measured variables 
has been the beginning point of the design of control and 
monitoring strategies. Usually, the variable whose estimate 
is desired is the one that is measured. Even the emerging 
field of parameter estimation for online optimization relies 
in practice on the information obtained through traditional 
monitoring. Analytical redundancy can provide better and 
more reliable estimates. We therefore call key variables 
those for which estimators are of interest for control, moni- 
toring, and parameter estimation. Thus, in general, mea- 
sured variables are not necessarily the same as the key 
variables. 

The concept of key variable is already used in inferen- 
tial control [60,61] and procedures for the selection of key 
variables in control are reviewed elsewhere (see [72] or [85]). 
Several techniques have been developed to establish the exis- 
tence of a fault (observability), determine its nature (diag- 
nosis), and take appropriate action (alarm and/or shutdown), 
while assuming that the instrumentation is free of biases 
and is well calibrated. A few books cover all this in detail 
[36,39,75]. These techniques range from fault trees, cause- 
effect digraphs, neural networks, knowledge -based systems 
to model-based statistical approaches like principal compo- 
nent analysis and PLS. The field is intertwined with the data 
reconciliation/rectification field to the extent that the same 
techniques are used sometimes in both fields under different 
names. A detailed account of all these issues can be found in 
the books by Bagajewicz [11,12]. 


INSTRUMENTATION DESIGN GOALS 

Several goals exist for every designer of a network of instru- 
ments in a plant: 

• Cost (including maintenance cost) 

• Estimability: a term coined to designate the ability 
to estimate a variable using hardware or software 
[20], The term is a generalization of the concepts 
of both observability and redundancy (section 
on Design of Sensor Networks for Process Fault 
Resolution provides more specific definitions of 
these terms) 

• Precision 

• Reliability and availability of the estimates of the key 
variables, which are function of the reliability/avail- 
ability of the measurements 

• Gross error robustness (see [10]): this encompasses 
three properties: 


- Gross error detectability: it is the ability to detect 
a certain number of gross errors larger than a 
specified size 

- Residual precision : it is the precision of estimates 
left when a certain number of instruments fail 

- Gross error resilience: it is defined as control over 
the smearing effect of undetected gross errors 

• Accuracy: it is the sum of precision and direct bias 
as well as software-induced bias obtained when data 
reconciliation is run 

• Ability to perform process fault diagnosis 

• Ability to distinguish sensor failure from process 
failure 

While all these goals seem equally important, cost has been 
the traditional objective function used in design, although 
other goals, such as precision, have also been used as objec- 
tive functions. The cost-based mathematical programming 
design procedure is the following: 

Minimize {Total Cost} 
s.t. 

Accuracy of Key Variables 

The total cost includes the maintenance cost, which regulates 
the availability of variables, a concept that substitutes reli- 
ability when the system is repairable. Earlier versions of this 
optimization model included precision, gross error robust- 
ness, and other goals separately. Latest work allows the use 
of accuracy as a replacement. In addition, there are two dif- 
ferent versions of accuracy: software accuracy or stochastic 
accuracy. 

The above design procedure does take into account the 
ability of a sensor network to detect process faults and logic 
for alarm systems. Even though important attempts are being 
made to address this issue [92], a model based on cost-effi- 
cient alarm design is yet to be produced. Likewise, the direct 
incorporation of control performance measures as additional 
constraints of a cost-optimal representation has been par- 
tially investigated, as we shall see below. Finally, methods 
for cost-optimal instrumentation design corresponding to the 
implementation of several other monitoring procedures, like 
principal component analysis (PCA), PLS, wavelet analy- 
sis, and neural networks, among others, have not been yet 
proposed. 

From the exclusive point of view of fault detection, the 
problem of the design of instrumentation is 

Minimize {Total Cost} 
s.t. 

• Desired Observability of Faults 

• Desired Level of Resolution of Faults 

• Desired level of Reliability of Fault Observation 

• Desired level of accuracy of the information gathered 
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The combination of both goals, that is, the design of a sensor 
network capable of performing estimation of key variables 
for production accounting, parameter estimation for online 
optimization as well as fault detection, diagnosis, and alarm 
(all of this in the context of detecting and assessing instru- 
ment biases and process leaks) together with attending the 
control needs using the same set of instruments is the long- 
term goal of this field. 

An alternative model is emerging. This model is based 
on the following mathematical programming design 
procedure: 

Maximize {Economic Value of Data - Total Cost] 

We notice that in this case there are no constraints. This 
is because accuracy is directly related to the Economic Value 
of Data. As we shall see below, the constraints that can be 
added to this model is to limit the number of sensors or the 
cost. The above procedure has been applied separately to 
three different perspectives, production accounting, fault 
detection, and identification and control. Only one paper by 
Bagajewicz et al. [16] made an attempt to study the design of 
one sensor network that will add the economic value from all 
perspectives. 

We now review design methods for the cost as well as the 
value paradigm. 


COST OPTIMAL ACCURATE SENSOR NETWORK 

Let c, be the cost of the candidate instrument to measure 
variable i, and let q he a vector of binary variables to denote 
whether a variable is measured (q ; = 1) or not (<:/, = ()). Then 
cost-optimal design of an accurate sensor network is obtained 
solving the following problem [10]: 



ieMi 


aj(q)<a* \/jeM P 
qi g {0, 1} VieM 


( 3 . 1 ) 


after data reconciliation and a* its threshold value. This 
accuracy can be obtained by adding the precision of the esti- 
mator after data reconciliation plus the maximum undetected 
induced bias in the key variable due to biases in the network, 
or through an average of possible outcomes obtained using 
Montecarlo simulations. Section on Data Reconciliation 
and Software Methods for Bias Detection describes accu- 
racy in more detail. We omit covering earlier models that 
target maximum precision in the network [1-3,11,12,53,54]. 
A mathematical connection between the maximum precision 
and the minimum cost representations of the problem exists 
[19]. More precisely, the solution of one problem is one solu- 
tion of the other and vice versa. 

The following small example, taken from Bagajewicz 
[10,1 1], which focuses on precision alone, is used to illustrate 
the procedure. 

EXAMPLE 

Consider the process flow diagram of Figure 3.1. The flow 
rates are F= (150.1,52.3,97.8,97.8). Assume that for each 
rate, flowmeters of precision 3%, 2%, and 1% are available 
at costs 800, 1500, and 2500 respectively, regardless of size. 
Precision is only required for variables F 1 and P 4 , that is, 
M P = {5 , 1 ,5 , 4 }, with o* = 1.5% and at = 2.0%. Two solutions 
are obtained featuring a cost of C= 3000. The corresponding 
meters are shown in Table 3.1. Assume now that the cost of 
the 3% sensors drops to 700. Then the optimal solution is no 
longer the one shown in Table 3.1. In this case, two solutions 
of equal cost (C=2900) are shown in Table 3.2. This solution 
is redundant and cheaper. 

Although precision is achieved using nonredundant net- 
works, gross errors are impossible to detect. Therefore, if at 
least one degree of redundancy is requested, that is, at least 
two ways of estimating each key variable, then there are two 


Si. 



where 

M P is the set of all variables for which accuracy con- 
straints are imposed 

M is the set of variables where measurements can be 
placed 

In the formulation of the objective function, it is assumed 
that there is only one potential measuring device with associ- 
ated cost c, for each variable (i.e., no hardware redundancy), 
but this condition can be relaxed [10,11], Furthermore, afq) 
represents the accuracy of the estimator of variable j obtained 


FIG. 3.1 

Two units— four streams example. 


TABLE 3.1 

Solutions of the Precision-Constrained 
Problem 


Solution Si 

S 2 (%) 

S 3 (%) 

s 4 (%) 

A — 

2 

2 

— 

B — 

2 

— 

2 


© 2012 by Bela Liptak 


3 Instrumentation in Processes and Automation 


75 


TABLE 3.2 

Solutions of the Precision-Constrained 
Problem (New Cost of Instruments) 

Solution S l (%) 

S 2 (%) S 3 (% ) S 4 (%) 

C 3 

D 3 

3 2 — 

3 — 2 


TABLE 3.3 

Solutions of the Precision-Constrained 
Problem ( Redundancy Required) 

Solution S x (%) 

S 2 (%) S 3 (%) S t (%) 

E 3 

F 3 

3 2 — 

3 — 2 


solutions with a cost of C=3100 (Table 3.3). This can be for- 
mally achieved adding a constraint requiring the estimability 
degree of each key variable to be two. The concept of degree 
of estimability is a generalization of the concepts of degrees 
of observability and redundancy [20] and it basically relates 
to the number of ways a variable can be estimated, regardless 
of whether it is measured or not. 

Bagajewicz [10] proposed a tree enumeration procedure 
as a means to solve the problem. In a recent work, it was 
shown that the problem can be in principle reduced to a con- 
vex mixed integer non-linear programming (MINLP) formu- 
lation [31]. Finally, Bagajewicz and Cabrera [14] showed that 
the model can be written as an mixed integer linear program- 
ming (MILP) problem, but also found some limitations to 
this approach. 

A series of techniques, which produce cost-optimal net- 
works subject to constraints on estimability, which were 
introduced by Bagajewicz and Sanchez [20], were extended 
to bilinear systems (typically those for which component 
and energy balances are added to overall mass balances) by 
Bagajewicz and Sanchez [21], 

PARAMETER ESTIMATION 

Considerable attention is being put nowadays to the issue 
of parameter estimation, especially in the context of the 
increasing popularity of online optimization. One of 
the early papers by Britt and Luecke [29] proposed the use 
of maximum likelihood principle to obtain parameters in 
implicit models. Reilly and Patino-Leal [80] initiated a line 
of work that bases parameter estimation on linearization but 
Kim et al. [42] proposed the use of nonlinear programming. 
The effect of data reconciliation and gross error detection 
in parameter estimation has been analyzed by MacDonald 
and Howat [52] as well as by Serth et al. [87] and Pages 
et al. [73]. Among the approaches based on linear alge- 
bra, Kretsovalis and Mah [44] proposed a combinatorial 


search based on the effect of the variance of measurements 
on the precision of reconciled values. Tjoa and Biegler 
[91] explored methods for the estimation of parameters in 
differential-algebraic equation systems. To select measure- 
ments for such goals, Krishnan et al. [45,46] presented a 
three-step strategy and Loeblein and Perkins [51] discussed 
the economics. 

The method proposed by Krishnan et al. [45] relies on 
a screening procedure that involves three steps: (1) A first 
step performs a structural analysis (singular value decom- 
position), which disregards measurements with little or no 
effect on the parameters, (2) a second step disregards mea- 
surements that have insignificant effect on the axis length of 
the confidence region of the parameter estimates, and (3) the 
last step determines the interaction between the parameter 
estimates by means of calculating a covariance matrix [24]. 
If the off-diagonal elements are too large, then the param- 
eters are highly interactive and therefore, any problem with 
the set of measurements that affects one parameter will 
also affect the other. The “best” set of measurements will 
have a small confidence region and lead to low interaction 
between the parameters. Unfortunately, this method does not 
take into account cost and does not offer a systematic proce- 
dure to make a final selection of the “best” set. In contrast, 
Bagajewicz [11] discusses linearization techniques in the 
context of cost-based minimization design procedures. 

PRECISION UPGRADE 

As pressure mounts to obtain more reliable and accurate esti- 
mates of variables and parameters, the use of data reconcili- 
ation techniques is the first response, but in most cases, the 
existing redundancy level is not enough to guarantee the level 
of accuracy required. Thus, one is faced with an upgrade 
problem. 

There are three possible ways of performing the upgrade 
of a sensor network: (1) by the addition of new instruments, 
(2) by the substitution of existing instruments by new ones, 
and (3) by relocation of existing instruments. 

Typically, addition of new instruments has been the 
response first considered. Kretsovalis and Mah [44] proposed 
a combinatorial strategy to incorporate measurements one at 
a time, to an observable system, but no constraints were con- 
sidered. Substitution and/or relocation are options that are 
sometimes cheaper. One example is the substitution of ther- 
mocouples by thermoresistances or their relocation. However, 
it is in the case of laboratory analysis where the options of 
substitution and/or relocation should be strongly considered. 
Since upgrading requires capital expenditure, it must be 
done on the basis of a cost-benefit analysis. The costs of the 
instrumentation are straightforward to obtain. However, the 
benefits need to be somehow quantified. In the case of data 
accuracy needs for accounting purposes, the benefit can be 
quantified as the decrease of lost revenue due to imprecise 
data. Describing it in simple terms, the larger the uncertainty 
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in the assessment of the amount of raw materials purchased 
and/or the products sold, the larger is the probability of lost 
revenue. Thus, every percent of accuracy of these flows can 
be assigned a monetary value. 

In the case of monitoring and parameter estimation for 
online optimization, an economic measure can also be devel- 
oped. First, in the case of monitoring, one can associate a 
revenue loss for product that is not in specification. Finally, 
in the case of online optimization, Loeblein and Perkins [51] 
discuss the use of measures of loss of economic performance 
due to offset from the correct optimum due to plant-model 
mismatches. In other cases, industry is looking less at the 
benefit in monetary terms and simply plots the increased pre- 
cision as a function of investment [3]. Although this approach 
is intuitive, its importance relies on the possibility of visual- 
izing the effect of instrumentation cost. 

The upgrading of a sensor network by the simple addition 
of instrumentation has to be done at minimum cost while 
reaching the goals of precision in key variables. Hence, the 
following minimum cost problem [22]: 


Min ^ ^c ik qit 

isMi ksKi 

S.t. 


Oj(q) < O* 

V; e M p 

J q* +N, < N* 

V; e M 

keKi 


ql e {0,1} 

V; e M,Vk e K { 


where 

Kf is the set of new sensors available to measure variable i 
Nj is the number of existing sensors measuring variable i 
c ik is the cost of the k-new sensor for measuring variable i 
N* is the maximum number of sensors that are allowed to 
be used to measure variable i 
q N is the vector of binary variables for new sensors 

The first constraint establishes a bound on precision, which is 
a function of the set of fixed existing sensors and new instru- 
mentation. The second constraint establishes an upper bound 
on the number of sensors used to measure each variable. This 
number is usually one for the case of flow rates, but it can be 
larger in the case of laboratory measurements of concentra- 
tions. In this case, these constraints can be lumped in one 
constraint for all concentration measurements to express the 
overall limitation of the laboratory on a daily basis, or frac- 
tion thereof. 

When maximum precision is requested, and cost is a con- 
straint, then a maximum precision problem can be used using 
a bound on the capital expenditure. The two problems are 
equivalent. 


EXAMPLE 

Consider the industrial heat exchanger network of Figure 3.2. 
It corresponds to a set of heat exchangers where crude is 
heated up using hot gas-oil coming from a column. In this 
case, the heat transfer coefficients for the heat exchangers 
are estimated in terms of temperature and flow rate measure- 
ments. The existing instrumentation is indicated in the figure 



F 1 = 217,019 lb/h r 9 = 266.9°F O Flowmeter 


Thermocouple 


Heat 

exchanger 

Area 

(ft 2 ) 

Ft 

C Ph (BT 
U/lb °F) 

Cp c (BTU/lb °F) 

U k 

500 

0.997 

0.6656 

0.5689 

u 2 

1100 

0.991 

0.6380 

0.5415 

U3 

700 

0.995 

0.6095 

0.52 


FIG. 3.2 

Industrial heat exchanger network. 
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TABLE 3.4 

Availability of New Instrumentation 



Flowmeters 


Temperature Sensors 


Standard 


Standard 


Stream 

Deviation (%) 

Cost 

Deviation (°F) 

Cost 

F 1 

3 

2250 

2/0.2 

500/1500 

f 2 

3 

2250 

2 

500 

F, 

3 

2250 

2 

500 

F, 

3 

2250 

2/0.2 

500/1500 

f 5 

3 

2250 

2 

500 

f 6 

3 

2250 

2 

500 

f 7 

3 

2250 

2 

500 

F s 

3 

2250 

2 

500 

f 9 

3 

2250 

2/0.2 

500/1500 


TABLE 3.5 

Maximum Number of 
Instruments for the Heat 
Exchanger Network 

Variable 

N * 

Variable 

N* 

C 

1 

p 

2 

f 2 

1 

T 2 

1 

F, 

1 

T, 

1 

f 4 

1 

T, 

2 

Fs 

1 

Ts 

1 

f 6 

0 

t 6 

1 

f 7 

1 

P 

1 

F« 

0 

T s 

0 

f 9 

1 

t 9 

2 


(flowmeters are of 3% precision and thermocouples have a 
precision of 2°F). 

The standard deviations of heat transfer coefficients of the 
heat exchangers calculated using the installed set of instru- 
ments are [12.27 2.96 3.06] (BTU/h ft 2 °F). To obtain these 
values, all redundant measurements have been used. In order 
to enhance the precision of the parameter, new instruments 
should be added. In this example, hardware redundancy is 
considered. Furthermore, different types of new instruments 
are available to measure some temperatures. Data for new 


instrumentation is presented in Table 3.4, where cost and 
standard deviation are shown. The maximum number of 
allowed instruments for measuring each variable is given in 
Table 3.5, where a zero value indicates there is a restriction 
for measuring the corresponding variable. 

Table 3.6 presents results for the upgrade problem using 
the minimum cost problem. When there are two possible 
instruments to measure a variable, the type of instrument 
is indicated between parentheses in the optimal solution 
set. Thus, for example, T 4 (1) indicates that the first instru- 
ment available to measure the temperature in stream ,S’ 4 is 
selected, in this case a thermocouple with a precision of 2°F 
is selected. The weights for the maximum precision problem 
are assumed equal to one. One can notice that as precision 
requirements are made more and more stringent, the cost and 
the number of instruments added increase. Case 3 has three 
alternative solutions. 

RESOURCE REALLOCATION 

As mentioned earlier, in many cases, measurements can be 
easily transferred at no cost from one stream to another. This 
is the case of concentration measurements that are performed 
in the laboratory. Pressure gauges and thermocouples can 
also be transferred from one place to another at a relatively 
small cost. However, flowmeters are probably an exception. 
Since one can consider that the resource reallocation does 
not involve cost, or eventually neglects it, the minimum 
cost problem would reduce to a set of mixed integer nonlin- 
ear algebraic inequalities from which one solution is to be 
extracted. However, even in the case where cost is not con- 
sidered, one would like to minimize the number of changes. 
In addition, these reallocation costs may overcome the simple 
addition of new instrumentation. Therefore, any reallocation 
and upgrade program should consider the trade-offs between 
all these decisions. A mathematical programming represen- 
tation has been developed [22] before, and now it will be 
illustrated next. 

EXAMPLE 

Consider the flash tank of Figure 3.3. The number of exist- 
ing instruments to measure each variable (Nj) and the 


TABLE 3.6 

Results for the Minimum Cost Problem 


Case 

* 

Cft/, 

Ou 2 

* 

<*t/ 3 

Of/, 

Ou, 

Gu 3 

Cost 

Optimal Set 

i 

4.0 

4.0 

4.0 

3.6160 

1.9681 

2.7112 

500 

T 6 ( 1) 

2 

3.5 

2.0 

2.5 

2.7746 

1.6892 

2.3833 

1500 

T 2 (1) T 4 (1) T 6 (1) 

3 

3.0 

1.5 

2.5 

2.7230 

1.4972 

2.2844 

6500 

F 2 F 3 T 2 (1)T 4 (1)T 6 (1)T 9 (1) 
F 2 F 4 T 2 (1) T 4 (l) T 6 (l) T 9 (l) 
F 3 F 4 T 2 (1) T 4 (l) T 6 (l) T 9 (l) 

4 

3.5 

2.0 

2.0 

— 

— 

— 

— 

— 
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TABLE 3.9 

Flash Drum: Results for Sensor Reallocation and Upgrade 

Case 

0 * 

a 

Cost 

Reallocations 

New 

Instruments 

i 

oo 

0.00438 

— 

— 

— 

2 

0.0038 

0.00352 

too 

— 

p 



0.00347 

too 

yi to y 2 

p 

3 

0.0033 

0.00329 

2800 

yi to y 2 

y 2 P 


FIG. 3.3 

A flash unit. 


corresponding standard deviations (q ik ) for different types of 
installed instruments are included in Table 3.7. In this exam- 
ple, the mass fractions of all the components of a stream are 
measured on line but a laboratory analysis may be done as a 
second alternative to know their values. 

The standard deviation of the vaporization efficiency 
coefficient estimated using the existing instrumentation is 
0.00434. As this value is not satisfactory, a reallocation and 
possibly an incorporation of new flowmeters and laboratory 
composition analysis are proposed. The sets of streams from 
which measurements can be transferred to other streams is 
M t ={ 1,2,3} and the set of new locations is M R = {1,3, 4, 6}. 
The costs of feasible reallocations between sets M T and M h . 
are given in Table 3.8. For example, the cost of relocating a 
sensor from h\ to F 2 is 80. Infinite costs are used for forbid- 
den reallocations. Typically, this relocation pattern is con- 
structed based on engineering judgment. 


TABLE 3.7 

Installed and New Instrumentation Data 


New Available 


Index 

Variable 

N, 

Installed o i k 

c i.k 


i 

F 

2 

2.5 

2.5 

350 

2 

2 

?! 

2 

0.015 

0.01 

2700 

0.01 

3 

f 2 

1 

1.515 

— 

350 

1.48 

4 

yi 

1 

0.01 

— 

2700 

0.01 

5 

F 

1 

1.418 

— 

400 

1.38 

6 

A 

1 

0.01 

— 

2700 

0.01 

7 

P 

1 

14 

— 

100 

14 


TABLE 3.8 

Flash Drum: Costs of 
Relocation 



F, 

F 

y 2 

y 3 

F 

OO 

80 

oo 

oo 

y t 

oo 

OO 

0 

50 

F 

80 

oo 

oo 

OO 


Table 3.9 shows the results. The first row represents the 
case for the existing instrumentation. A reduction of the stan- 
dard deviation from 0.00438 to 0.00347 results if the labora- 
tory analysis for the feed stream is relocated to the liquid 
stream and a pressure sensor is added. The cost of this case 
is 100. Higher precision is obtained by the reallocation and 
addition of instruments. For a* = 0.0031, no reallocation and 
instrument addition can meet this goal. 


RELIABLE SENSOR NETWORKS 


Bagajewicz [11,21] discusses the concepts of availability and 
reliability for the estimates of key variables and distinguishes 
them from the same concepts applied to instruments. The 
former are called estimation reliability/availability whereas 
the latter are called service reliability/availability. Ali and 
Narasimhan [4] proposed to use the system reliability given 
by the minimum estimation reliability throughout the whole 
network as the objective function. In other words, the reliabil- 
ity of the system is defined by its weakest element. To address 
the fact that their representation does not control cost, they 
proposed to limit the number of sensors to the minimum pos- 
sible that will still guarantee observability. Although their 
procedure does not guarantee global optimality, it produces 
good results. Bilinear networks are also discussed by Ali and 
Narasimhan [5] in detail. Finally, genetic algorithms are suc- 
cessfully used by Sen et al. [86] not only for reliable networks 
but also for a variety of other objective functions. 

The cost-based representation for the design of the sensor 
network subject to reliability constraints is [21] 


Min 


X Ci q ' 

V/eMj 


R l k (q)>Rt VkeM R 

g,e{0,l] VieM 


(3.3) 


where 

Rl(q) is the reliability of variable k 
R * is the threshold 

M r is the set of variables whose reliability is to be 
constrained 
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The reliability of each variable is calculated using the fail- 
ure probabilities of all the sensors participating in the corre- 
sponding cut sets (material balances). If all sensors have the 
same cost c and that N is the number of sensors, one obtains a 
problem where the number of sensors is minimized. Finally, 
the representation due to Ali and Narasimhan [4] can be put 
in the form of a minimum cost problem. The details of such 
equivalency can be found in the article by Bagajewicz and 
Sanchez [21], where examples are shown. 


REPAIRABLE SENSOR NETWORKS 

When repairs are not present, the service availability of a 
sensor is equal to its service reliability. In addition, the fail- 
ure rate has been considered in a simplified way as a con- 
stant. However, in the presence of repairs, failure is no longer 
an event that depends on how many hours the sensor survived 
from the time it has been put in service. It is also conditioned 
by the fact that due to preventive or corrective maintenance, 
the sensor has been repaired at a certain time after being put 
in service. These events condition the failure rate. We thus 
distinguish unconditional from conditional events in failure 
and repair. These concepts are important because sensor 
maintenance cost accounts for nearly 20% of all mainte- 
nance cost [57]. Its reduction or containment is therefore 
essential. The connection between failure rate, repair rate, 
and the expected number of repairs as well as illustrations 
of the impact of maintenance on sensor network design are 
described by Bagajewicz [11]. 


EXAMPLE 

We now show a design of a sensor network for the simplified 
ammonia network (Figure 3.4). In this example, sensors for 
each stream may be selected from a set of three instruments 
with different precision, purchase cost, and failure rate. 
These data are included in Table 3.10. 

Maintenance corrective costs are evaluated considering 
spare part cost and manpower cost of 10 and 40, respec- 
tively, a life cycle of 5 years, and an annual interest of 6%. 



FIG. 3.4 

Simplified ammonia plant network. 


TABLE 3.10 

Instrumentation Data 


#i 

#2 

#2 

Purchase cost 

350 

250 

200 

Precision (%) 

1.5 

2.5 

3 

Failure rate (failure/year) 

0.3 

0.6 

0.7 


TABLE 3.11 

Constraints of the Optimization Problem 


Stream 

Precision 
Requirements (%) 

Residual Precision 
Requirements (%) 

Availability 

Requirements 

Ft 

— 

— 

0.9 

f 2 

1.5 

2.0 


f 5 

2.5 

3.0 


f 7 

— 

— 

0.9 


Constraints of precision, residual precision, and availability 
are included for only two flow rates. The limits on these vari- 
ables are presented in Table 3.11. The repair rate of instru- 
ments, which is a parameter as a characteristic of the plant in 
consideration, has been varied between 1 and 20. The results 
of the optimization problem are presented for each case in 
Table 3.12. 

In the first case, the repair rate is comparatively low. 
Consequently, the availability of instruments in the life 
cycle is also relatively low. To satisfy the availability of key 
variables, the optimal solution includes a set of six instru- 
ments. Three of these instruments are of type 1, sensors 
of low failure rate, high precision, and high cost. For this 
reason, precision and residual precision constraints are not 
binding. 

When the repair rate is 2, an optimal solution exists, 
which consists of five instruments. Two of these instruments 
are of type 1 and the rest are of type 3. Consequently, the 
total instrumentation cost decreases. 

A lower instrumentation cost is obtained for a repair 
rate equal to 4. Even though sensors are located on the same 
streams as in the previous case, one sensor of higher failure 
rate is installed on stream 5 8 . This occurs because the repair 
rate is now higher, making the constraint on availability of 
variable S-, not binding. 

The results of the last case show that the influence of 
availability constraints decreases for high repair rates. The 
cost increases because of the effect of increasing the repair 
rate (from 4 to 20) in the maintenance cost model. 

As a conclusion, the repair rate has a direct influence on 
the availability of a variable. If the repair rate is high, the 
design follows the requirements of precision and residual 
precision constraints. Thus, the availability of a variable may 
be a binding constraint for lower repair rates. In this situa- 
tion, cost may increase because it is necessary to incorporate 
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TABLE 3.12 

Optimization Results for the Simplified Ammonia Process Flowsheet 


Repair Rate 

Measured Variables 

Instrument 
Precision (%) 

Cost 

Precision (%) 
(F 2 ) 

(F 5 ) 

Availability 

(Ff 

(F 7 ) 

i 

F\ F t F 5 F 6 F 7 F g 

311132 

2040.2 

0.8067 

1.2893 

0.9021 

0.9021 

2 

Ft F 5 F 6 F 7 F g 

33 13 1 

1699.8 

0.9283 

1.9928 

0.9222 

0.9062 

4 

F i F 5 F 6 F 7 F s 

3 3 13 3 

1683.7 

1.2313 

1.9963 

0.9636 

0.9511 

20 

Ft F 5 F 6 F 7 F g 

3 3 13 3 

1775.2 

1.2313 

1.9963 

0.9983 

0.9969 


more instruments to calculate the variable by alternative 
ways. 

ROBUST SENSOR NETWORKS 

When data reconciliation is performed, systematic errors 
(biases) produce a smearing effect in all the data so that the 
whole data set has to be discarded. Without analytical redun- 
dancy, no smearing exists, but then there is no other way than 
direct inspection to detect sensor failure. 

A robust sensor network provides meaningful values 
of precision, residual precision, variable availability, error 
detectability, and resilience. These five properties encompass 
all the most desired features of a network. Indeed, precision 
and residual precision guarantee that data are always of the 
desired quality. Residual precision controls redundancy and 
reliability and allows the data to be most of the time at hand. 
Finally, error detectability and resilience make sure that the 
probability of data being free of gross errors is large. While 
the two first properties are deterministic and depend directly 
on the quality of the instrument selected, the other three are 
statistical in nature. It is expected that more properties will 
be added to define robustness in the future. For example, for 
neural networks, wavelet analysis, PCA, PLS, and other tech- 
niques for process monitoring, very little analysis regarding 
goals that a sensor network design or upgrade procedure can 
undertake was investigated. It is expected that in the next few 
years, such methods will emerge. Robustness can be incor- 
porated explicitly in the form of constraints to the minimum 
cost problem as discussed before [11]. We illustrate them 
next. 


EXAMPLE 

We now add residual precision capabilities to the example 
of Figure 3.1. We formally define residual precision of order 
k for stream i, \\ij(k), the precision that the estimate of the 
flow rate of stream i will have if any k measurements of the 


network are eliminated. Consider now that residual precision 
of order k= 1 is added to variables .S', and ,S' 4 as follows: 
\|/f (1) = 1.5% and \]/*(l) = 3%. The solution is to put sensors 
of precision 2%, 3%, 3%, 3% in S) through S 4 , respectively. 
The cost is C= 3900. Assume now that residual precision is 
requested to the same level as precision. Then two alternative 
solutions with cost C=5500 are obtained (Table 3.13). 

Not only the cost is higher but also there is one more 
degree of redundancy. For larger problems, the number of 
alternatives will increase, requiring new criteria to further 
screen alternatives. 

We now turn to adding error detectability to this exam- 
ple. As the error detectability increases, so does the precision 
of the sensor network. However, it can also make the design 
infeasible. Consider adding an error detectability level of 
k ; , = 3.9 to the problem. Requiring such a value of k d means 
that the network should be able to detect gross errors of 
size 3.9 times the standard deviation of the measurement or 
larger in any measured variable. A statistical power of 50% 
is assumed in the algorithm used to detect the gross errors. 
Two solutions from a set of only 4 feasible solutions are found 
with cost C=4800 (Table 3.14). 


TABLE 3.13 

Solutions of the Residual 
Precision-Constrained Problem 

S x (%) S 2 (%) S,(%) S 4 (%) 

1 2 2 — 

1 2 — 2 


TABLE 3.14 

Effect of Error Detectability 
Constraints (k d = 3.9 ) 

S,(%) S 2 (%) S,(%) S t (%) 

1 3 — 2 

1 3 2 — 
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If an error detectability of k d = 3.4 is requested, the problem 
has only one solution, namely, sensors with precision of 1%, 
3%, 1%, 1% in .S' | through S 4 , respectively, with a cost C=8300. 

A resilience requirement means that the network will be 
able to limit the smearing effect of gross errors of size a cer- 
tain size without affecting the other properties. If resilience at 
a level of three times the standard deviation for all measure- 
ments is added, then the solution is to put sensors with preci- 
sion 1%, 3%, 1%, and 1% in 5) through S 4 , respectively, with 
a cost of C= 8300. Relaxing (increasing) the resilience levels 
while maintaining the error detectability at the same level may 
actually lead to solutions of higher cost, even to infeasibility. 

Example of Accuracy-Constrained Networks 

Consider the process flow sheet is shown in Figure 3.5, the 
data for the example is shown in Table 3.15. Sensor failures 
are assumed to occur with probability 0.2 (third column), 
the biases associated to these failures are assumed to follow 


normal distributions with zero means (fourth column) and 
standard deviations given in fifth column. For this example, 
the results are shown in Table 3.16 for various threshold val- 
ues and key variables. 

The desired value of accuracy decreases (from 4.0 to 1.5), 
which requires more sensors to be used. In design case a, only 
observability is required for the two key variables S l and S 5 . 
In design cases b and c, the obtained optimal solutions render 
redundancy of level one for key variables (the two key vari- 
ables are still observable if one removes any one sensor out of 
the three-sensor solutions). When a smaller accuracy thresh- 
old is required (design case d), both redundancy (of key vari- 
ables) and gross error detection capability of the network are 
required; hence more sensors need to be used. It can be seen 
that the optimal solution in this design is the same as the solu- 
tion obtained in design c where both estimation redundancy 
and gross error detection capability are required. Design case d 
is the extreme case where the required accuracy threshold is so 
small that all sensors need to be used to meet the requirement. 



FIG. 3.5 

Flowsheet to illustrate an accuracy-constrained design. 


TABLE 3.15 

Data for the Example in 

Figure 3.5 



Stream 

Flow 

Rates 

Costs 

Probability 
of Failure 
of Sensor 

Mean of 
pdf of Bias 

STD of pdf 
of Bias 

S , 

2 

55 

0.2 

0 

8 

s 2 

2 

40 

0.2 

0 

11.2 

s 3 

2 

60 

0.2 

0 

11.2 

s 4 

2 

50 

0.2 

0 

1.6 

s 5 

2 

45 

0.2 

0 

9.6 

s 6 

2 

55 

0.2 

0 

1.6 

s 7 

2 

60 

0.2 

0 

8 


DESIGN OF ECONOMIC VALUE-OPTIMAL 
SENSOR NETWORKS 

As stated above, the sensor network design formulation is the 
following: 

Max{V( 4 r)-c(tf)} 

s.t. 



i= 1 


c {(( ) 5s c max J 

where 

M, is the number of candidate sensors (number of process 
variables under considerations) 

V(q ) is value of sensor network (function of measurement 
locations q ) 
c(q ) is cost of sensors 
/V milx is limit on number of sensors 
c max is the limit on the budget 


TABLE 3.16 

Results for the Example 





Case Study a 

b 

c 

d 

e 


Key variables 

S) and S 5 

S) and S 5 

Si and S 5 

Si and S 5 

Si and S 5 

Accuracy 

thresholds 

4 

3 

2 

1.8 

1.5 

Accuracy value 

a SI = 3.36 

a sl = 2.22 

«si = 1-90 

a SI = 1.65 

«si = 1-499 


a S5 =2.85 

a s5 = 1-99 

o S5 = 1.81 

« S5 = 1-49 

a s5 = 1-27 

Measured variables 

s»s 6 

Si, s 5 , s 6 

Si, S 6 , S 7 

Si, S 5 , S 6 , S 7 

All variables 

Cost 

110 

155 

170 

215 

365 
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For the production-accounting perspective, the value is given 
by the difference between the downside expected financial 
loss when no sensors are used (a reference value) and the 
downside expected financial loss corresponding to the sensor 
network: 


V(q) = {DEFL ( no sensor) - DEFL ( with sensors)} (3.5) 

In turn, the downside expected financial loss is obtained as 
follows: Because of the random errors of the estimator of a 
product stream flow rate, there is a finite probability that the 
estimator is above the target but in fact the real flow is below 
it. In such situation and under the assumption that the opera- 
tors do not make any correction to the production through- 
put set point when the estimator suggested that the targeted 
production has been met or surpassed, the production output 
will be below the target and financial loss occurs. The finan- 
cial loss under simplified assumptions of negligible process 
variations and normal distributions of the process variation 
and the measurements was found to be DEFL = 0.19947* T* 

E K sp a p where K is the cost of the product (or the cost 

p 

of inventory) for flow p, T is the time window of analysis [17] 
and the precision of the estimator a p . 

Using the same concept of downside financial loss, 
Bagajewicz [13] extended the theory of economic value of 
precision to include the effect of (induced) bias, namely, the 
economic value of accuracy. The expression for financial loss 
DEFL considering bias is given by [13] 


DUEL = ^ 0 DEFL° + ^ ^DEFL 1 1 + ^ '¥j u2 DEFL 2 ^ 

i il,i2 

+ - Y V^.jnDEFL 11 I 


(3.6) 


In this equation, V F" 1 , 2 iW and DEFL N \ are the average 
fraction of time the system is in the state containing n gross 
errors q, i 2 ,...,i N and its associated financial losses, respec- 
tively. Detailed expressions and procedures to calculate the 
financial loss for system containing n biases q, i 2 ,. ,.,i N can be 
found in Nguyen Thanh et al. [70]. 


EXAMPLE 

We use the same network discussed in Figure 3.5 and data 
of Table 3.15, and there is no addition of adverse conditions 
on the number of sensors used or on the cost. The time win- 
dow of analysis T is 30 days, the cost of product K s (or cost 
of inventory) for the two key variables .S', and S 5 is shown 
in Table 3.17. The financial losses of the optimal sensor net- 
works are shown in row 6 of the table 


TABLE 3.17 

Results for Value-Based Design 

Case Study 

a 

b 

c 

d 

Key variables 

S, and S 5 

Sj and S 5 

Sj and S 5 

S 7 snd S 5 

K s value 

^ t ,=2 

K S1 = 10 

K, i = 30 

o 

kO 

II 

fc<f 


K s 5=2 

K, 5 = 10 

o 

<N 

II 

II 

CD 

o 

Measured variables 

s„s 6 

Sn S 6 , S 7 

*^5’ *^7 

All 

Sensors cost 

110 

170 

215 

365 

Financial loss 

78.6 

219.9 

451.8 

824.9 


When the cost of product K s increases, financial loss 
increases. If K s is small, the cost factor dominates financial 
loss factor and the optimal network contains few sensors. In 
the opposite site, if K s is large, the financial loss term domi- 
nates the cost term and the optimal network contains large 
fraction of candidate sensors so as to minimize financial loss. 

In design case a, K s is small, cost needs to be minimized, 
and the optimal network contains only enough sensors to 
guarantee observability of key variables. 

In design case d, K s is large, financial loss needs to be 
minimized, and the optimal network contains all sensors. 
This is an extreme case. 

In design cases b and c, K s is moderate, the optimal net- 
works contain enough sensors, which can guarantee some 
degree of estimation redundancy and gross error detection 
capability. The optimal networks in these two design cases 
are the same as the networks obtained in the two design 
cases c and d, which have good process monitoring capa- 
bility (good accuracy value and good gross error detection 
capability). 


DESIGN OF SENSOR NETWORKS FOR 
PROCESS FAULT DIAGNOSIS 

Whilst the basic goal of monitoring systems is to provide a 
good estimate of the state of the system, alarm systems are 
design to alert about process malfunction. In turn, process 
faults, which typically are rooted in some unit, propagate 
throughout the process, altering the readings of instruments 
(pressures, temperatures, flow rates, etc.). Thus, these sensors 
should be able to determine departures from normal opera- 
tion. In this sense, this task is different from that of gross 
error detection, which concentrates on instrument malfunc- 
tion. As a consequence, the discrimination between instru- 
ment malfunction and process fault is an additional task of 
the alarm system. As a consequence, the problem of design- 
ing an alarm system consists of determining the cost-optimal 
position of sensors, such that all process faults, single or mul- 
tiple and simultaneous, can be detected and distinguished 
from instrument malfunction (biases). 

A process fault is a departure from an acceptable range 
of operation or “degradation from normal operation condi- 
tions, and includes symptoms of a physical change (such as 
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deviations in measured temperature or pressure) as well as 
the physical changes themselves (scaling, tube plugging, etc.) 
and deviations in parameters such as heat transfer coefficients 
[96].” Faults originate in a process and propagate to a set of 
sensors. These sensors are also subject to faults themselves. 
These faults on the sensors are either biases or catastrophic 
failures. In the latter case, the fault detection is compromised, 
whereas in the former case, the process fault can go either 
undetected or false alarms may be induced. Therefore, a good 
alarm system should be able to filter the disturbances affecting 
the sensors as well as the gross errors induced by their faults. 
The next step is the process fault detection itself, although 
these two steps can be performed simultaneously. This pro- 
cedure needs also to have some capabilities of distinguishing 
process disturbances from real faults. Once a process fault 
has been identified, the final step consists of taking correc- 
tive actions, or determining a shutdown of the process. This 
is performed by implementing an alarm logic. Fault detection 
and diagnosis has been addressed in several books [36,39,75]. 
Bagajewicz [11] offers a review of the different approaches. 

The first attempt to present a technique to locate sensors 
was done by Lambert [49], where fault trees are used based 
on failure probabilities. Failure probabilities are hard to 
assess and fault trees cannot handle cycles and the construc- 
tion of the tree is cumbersome for large-scale systems. Due 
to these limitations, the technique has not been developed 
further. Raghuraj et al. [79] proposed an algorithm for one 
and multiple fault observability and resolution. Their design 
procedure is not based on cost and does not include other 
features, like filtering of sensor failures. Other aspects that 
have not been considered yet are failure probabilities, sensor 
failure probability, and severity of particular faults and the 
cost of the sensors. 

Raghuraj et al. [79] used directed graphs (DG), that is, 
graphs without signs. The arcs of the DG represent a “will 
cause” relationship, that is, an arc from node A to node B 
implies that A is a sufficient condition for B, which in general 
is not true for a Signed DG, where an arc represents a “can 
cause” relationship. The strategy used to solve the problem is 
based on identifying directed paths from root nodes where 
faults can occur to nodes where effects can be measured, 
called the observability set. Of all these paths, the objective is 
to choose the minimal subset of sensors from the observabil- 
ity set that would have at least one directed path from every 
root node. Raghuraj et al. [79] proposed a greedy search, 
which is later modified to remove redundant members of the 
observability set by a “backtracking algorithm.” They do 
not guarantee optimality, that is, a minimum set may not be 
found by this algorithm. It is likely that a formal mathemati- 
cal programming version of the problem will accomplish it. 

Consider the DG of Figure 3.6; all the nodes of the 
observability set are in the top row and the root nodes are 
in the bottom row. The problem of choosing the minimum 
number of sensors (nodes in the top row) that would cover 
all the root nodes is the well-known “minimum set covering 
problem [74].” All the root nodes are said to be “covered” if 



FIG. 3.6 

Bipartite Graph. 

a directed path exists from every root node to at least one of 
the nodes of the observability set. 

The solution to this problem obtained by Raghuraj et al. 
[79] is that nodes C 7 and C 8 constitute the minimal set. Indeed, 
R l and R 2 are observable from C 7 and the rest from C 8 . 

DESIGN OF SENSOR NETWORKS FOR 
PROCESS FAULT RESOLUTION 

For a fault-monitoring system to be useful in practice, it 
should not only be able to observe all the faults but also 
identify their location and number to the maximum extent 
possible. This is called maximum fault resolution. Therefore, 
a sensor network for maximum fault resolution is such that 
each fault has one and only one set of nodes from which it 
is observable. 

Consider three fault nodes R t , R 2 , and R t (Figure 3.7). 
Clearly, if only one fault is expected to occur at a time, then 
the set [Cj, C 3 ] would be adequate to distinguish between the 
three faults. Indeed, a fault in If is reflected in C,, but not in C 3 . 
Similarly, a fault in R 2 is reflected in C 3 , but not in Cj. Finally, 
a fault in R } is reflected in both Cj and C 3 simultaneously. 

Raghuraj et al. [79] proposed a procedure to obtain such 
sets of nodes for maximum resolution. They also show that 
the sensor-location problem for multiple faults and maximum 
resolution can be solved as an extension of the single -fault 
assumption problem. The application of these algorithms to 
a continually stir tank reactor (CSTR) and a fluid catalytic 
cracking (FCC) unit are presented by Raghuraj et al. [79]. 

Bhushan and Rengaswami [25,26] made considerable 
progress in the study of several important aspects of the prob- 
lem. Later, Bagajewicz et al. [15] presented a methodology to 
obtain cost-optimal networks for fault observation and mul- 
tiple fault resolution. Moreover, Bhushan and Rengaswamy 
[27,28] presented a methodology for fault detection, which 
minimizes cost and takes into account the reliability of the 
network. 



Maximum fault resolution. 
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DESIGN OF SENSOR NETWORKS FOR PROCESS 
FAULT DIAGNOSIS— VALUE PERSPECTIVE 

Despite its importance, not too many attempts have been 
made to assess the value obtained from the ability to detect 
and observe faults. One such attempt is based on comput- 
ing the loss of profit that stems from non detecting/resolving 
different parametric faults [68]. The idea is quite general, 
however, the actual determination of the value function 
required some assumptions: exclude structural faults (equip- 
ment severe malfunction) and include only parametric faults 
(deviation from normal operation), single fault occurrence, 
that is not more than one fault at a time, instant correction 
of detected faults. They also used a probabilistic descrip- 
tion of faults and profit that is a certain function over certain 
period of time T of the deviation from normal values. Finally, 
they defined value as the expected profit. Because the profit 
function is monotone nonincreasing, value is assessed as the 
upper bound of profit and computed only in the first period. 
They illustrate the approach using a CSTR example. 

An alternative approach is to resort to Montecarlo simula- 
tions, similar to those performed by Nguyen et al. [71] to assess 
the effect of preventive maintenance frequency in the overall 
plant economic losses due to faults. Bagajewicz [12] also dis- 
cuss the optimization of the maintenance-related manipulated 
variables (preventive maintenance frequency, prioritization of 
work policies, inventory levels, and available labor). 

DESIGN OF SENSOR NETWORKS FOR AUTOMATION 

When one considers the broader perspective of process auto- 
mation, the data collected by a sensor network will be used to 
make decisions about plant operation. This structure of infor- 
mation feedback and process response suggests an alternative 
definition of sensor network performance from the automation 
point of view. Specifically, performance of the sensor network 
should be measured with respect to its ability to improve the 
functionality of the controller. That is, rather than specify 
the desired level of data collection accuracy, one should spec- 
ify the desired closed-loop control system performance. One 
measure of control system performance is the variability of key 
process variables. For example, the standard deviation of a con- 
trol variable (CV ) in a control system, around its set point could 
be a measure of closed-loop variability. Similarly, the standard 
deviation of a manipulated variable (MV) could also be used to 
indicate closed-loop performance. In many cases, MV variabil- 
ity will be more important than the CV criteria, as the MV cri- 
teria will often reflect physical limits of the actuator hardware. 

An additional issue with the automation perspective is 
the selection of the feedback element (or controller). In the 
case when only data collection is performed (e.g., for moni- 
toring or production accounting purposes), the construction 
of a data processing element (for conversion of data into 
information) is based solely upon the physical parameters of 
the process and selected instrumentation. From there, one 


can simply compare the quality of information obtained with 
the desired performance. However, in the process automa- 
tion case, there are additional degrees of freedom in relation 
to how this information will be used to influence plant opera- 
tion. Consider the simple example of a proportional-integral 
(PI) controller in a single-input, single-output servo-loop. If 
the sensor used to measure the CV is very noisy, then one 
would likely detune (or reduce the gain of) the controller to 
reduce the impact of this noise on the MV. However, if the 
CV sensor is of high quality, one would likely specify an 
aggressive tuning to the controller, to make the best use of 
the collected information. In general, the relation between 
data collection quality and controller tuning parameters is a 
feasibility question parameterized by the variability criteria 
discussed in the previous paragraph. That is, for a given set 
of variability criteria and a candidate sensor network, there 
will be a family of controllers capable of meeting the criteria. 

The topic of measurement device selection for control 
system applications has a long and distinguished history. The 
subject was first addressed by Kushner [48], Meier et al. [58], 
and Athans [9], where focus was on measurement time selec- 
tion for communication systems. Johnson [41] and Muller and 
Weber [65] both focused on the subject of selecting a sensor 
array to maximize the observability of a dynamic process. 
Yu and Seinfeld [98] proposed a scheme for the selection of 
measurement locations within a dynamic distributed param- 
eter system using objective of minimizing estimation error 
variance. (Similar approaches can be found in Harris et al. 
[37] as well as in the citations of the survey paper by Kubrusly 
and Malebranche [47]). Other early efforts from the chemi- 
cal process control literature include Mellefont and Sargent 
[59], Morari and O’Dowd [62], Morari and Stephanopoulos 
[63,64], and Romagnoli et al. [84], In the context of control 
structure (CV and MV) selection, an enormous body of liter- 
ature has been created; see, for example, Morari et al. [63,64], 
Hovd and Skogestad [40], Narraway and Perkins [69], Lee et 
al. [50], Loeblein and Perkins [51], Cao and Rossiter [30], 
Wisnewski and Doyle [97], Heath et al. [38], and Van de Wal 
and de Jager [93]. Recent contributions not in the context 
of control structure selection include; Van de Wouwer et al. 
[94], Faulds and King [35], Antoniades and Christofides [7], 
Chmielewski et al. [31], Muske and Georgakis [66], Alonso 
et al. [6], Peng and Chmielewski [76,77], Singh and Hahn 
[88], and Armaou and Demetriou [8]. 

Based on the above discussion, an automation-based, 
cost-optimal sensor network can be obtained by solving the 
following problem [77]: 
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where 

q is the vector of binary decision variables indicating the 
placement of a sensor at location i ( M the set of location 
where measurements can be placed) 

L is a matrix of tuning variables for the feedback controller 

Furthermore, oy’(q,L ) represents the variability of the per- 
formance variable z,- (as a function of the candidate sensor 
network and controller tuning values, q and L) and of’ * * is 
the desired bound (or criterion) on the variability of z,-. Details 
on the functionality and computational aspects of Cj(q,L) 
can be found in Peng and Chmielewski [77]. 

EXAMPLE OF AUTOMATION-BASED COST-OPTIMAL 
SENSOR NETWORKS 

Consider the process depicted in Figure 3.8. Assume there are 
two possible locations one could place temperature sensors; 
the furnace exit and the reactor exit. Additionally, assume 
the availability of three sensor types, each with a different 
precision and cost. Type 1 has a precision of 0.175 K and an 
annualized cost of $250, type 2 has precision of 0.1 K and a 
cost of $450, and type 3 has precision of 0.075 K and a cost 
of $800. The set of performance variables (z.-'s) is selected 
to be all of the state and manipulated variables indicated in 
Figure 3.9. 

The dashed lines in the plots of Figure 3.9 indicate the 
desired bounds with respect to variability for each perfor- 
mance variable (distance between this dashed line and the 
operating point, center dot, is the allowed standard deviation 
for that variable). The ellipse extending beyond the dashed 
box of the top left plot of Figure 3.8 is the expected dynamic 
operating region (EDOR) under open-loop conditions 
(assuming no feedback). The fact that this ellipse extends 
beyond the desired performance bound (the dashed line) 


indicates that the open-loop configuration is not feasible. The 
other ellipses of the figure represent the EDOR of a closed- 
loop system using a $450 sensor at the furnace exit. Since 
this configuration is also the solution to the design problem of 
Equation 3.7, the following statements are also true. Use of a 
less costly configuration (such as just one $250 sensor) would 
be infeasible in the sense that there does not exist a controller 
L such that the EDOR is contained within the performance 
bounds. In addition, while a more accurate configuration, 
such as the type 3 sensors at the furnace exit, will allow for 
a smaller EDOR, such a reduction will provide no benefit in 
the optimal cost formulation of Equation 3.7. These notions 
are illustrated in Figure 3.10. 

Within the process automation perspective, the notion of 
instrumentation network design can be expanded to include the 
selection and placement of actuators. Conceptually, such a for- 
mulation would be identical to that of Figure 3.8, with the only 
change being that the vector q includes binary variables indi- 
cating the attendance of actuators at various locations, in addi- 
tion to the sensors of the original formulation. Details of this 
formulation can also be found in Peng and Chmielewski [77]. 

The automation scenario also allows for an economic 
value-based formulation as explained before. In contrast to 
the cost-optimal formulation, this perspective defines eco- 
nomic value as a continuous function of EDOR size and 
shape. The basic idea being that a higher quality sensor (and/ 
or actuator) network will result in a smaller EDOR (or advan- 
tageously shaped EDOR) and thus higher value. Combining 
this value relation with a hardware cost relation, one can eas- 
ily define a design problem that is conceptually identical to 
that outlined here. Details on the value formulation for auto- 
mated process can be found in Peng and Chmielewski [76] as 
well as Bagajewicz [12], 

CONCLUSIONS 


CO, 0 2 out 
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Vent 

position 


PT 0 



Fuel feed rate 


Disturbance input: 

• Feed temperature 


Manipulated variables: 

• Reactant feed rate 

• Fuel feed rate 

• Vent position 



T F 






State variables: 

• Reactor temperature 

• Furnace temperature 

• Furnace 0 2 

• Furnace CO 


We presented an overview of cost-optimal and economic 
value-based formulations of the instrumentation network 
design and upgrade problem. The design of networks with 
constraints such as accuracy in key variables as well as 
other requirements, including its upgrade was illustrated. 
Additional features as reliability and the ability to handle 
instrument malfunction were also discussed. The design of 
instrumentation networks for fault detection and resolution 
was also briefly discussed. Finally, the very important prob- 
lem of sensor positioning for automation was discussed. An 
integrated perspective, outlined by Bagajewicz et al. [23], 
where a single set of sensors is used for all monitoring, fault, 
and automation purposes, is yet to be developed. 


Nomenclature 


FIG. 3.8 

Furnace-reactor process. 


a Vector of accuracy values 

a* Threshold of desired accuracies 


© 2012 by Bela Liptak 


86 Process Control and Automation 



8 8.5 9 9.5 10 10.5 11 11.5 12 

Fuel feed rate (bbl/day) 


0.09 0.095 0.1 0.105 0.11 

Vent position (%) 


Cost of product (or the cost of inventory) for flow p 
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Pressure 

Binary vector indicating whether sensors are 
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Binary vector for new sensors 

Reliability of variables 

Threshold reliability 

Stream i 

Temperature 
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Value of sensor network define by the vector q 

Mole fraction 

Precision of estimates 

Precision threshold 

Residual precision vector of order k 

Fraction of time the system has k simultaneous 
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